echo on  
%  calcolo delle seguenti funzioni equivalenti
% 1) x*(sqrt(x^2+1)-x
% 2) x*sqrt(x^2+1)-x^2
% 3) x/(sqrt(x^2+1)+x)
pause
% facciamo il grafico delle funzioni
% nell'intervallo [10^0,10^9]
% per fare il grafico fissiamo alcuni
% punti sull'asse delle x e calcoliamo
% le tre funzioni in questi punti
pause
% con l'istruzione logspace calcoliamo
% 500 punti in scala logaritmica 
x=logspace(0,9,500);
pause
% calcolo le funzioni
% su tutto il vettore delle x
r1=x.*(sqrt(x.^2+1)-x);
r2=x.*sqrt(x.^2+1)-x.^2;
r3=x./(sqrt(x.^2+1)+x);
pause
%faccio il grafico in scala logaritmica
% sull'asse x
semilogx(x,r1,'r.',x,r2,'b.',x,r3,'g-');

